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Abstract. We report deep, high spectral resolution WSRT H I 21cm observations of four high latitude compact radio sources, 
that have revealed a new population of tiny, discrete clouds in the diffuse ISM, with peak optical depths r ~ 0. 1 -2%, H I column 
densities of 0.4-8 xlO 18 cirr 2 and core temperatures of 20-80 K. Imaging detections confirm these low column densities and 
imply linear core dimensions of a few thousand AU, assuming a distance of 100 pc. The physical origin of these tiny Hi 
structures and their distribution in the ISM is at present unknown. Further observations will be required to determine whether 
they are a ubiquitous component of the ISM. 
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1. Introduction 

Despite more than three decades of Hi 21cm absorption & 
emission studies, we still do not have a good understanding 
of the physical characteristics of the atomic component, one of 
the most important phases of the interstellar medium (ISM). 
Crucial observational lacunae include the morphology of neu- 
tral hydrogen "clouds", the time-scales of their formation and 
destruction (i.e. whether they are transient or relatively stable 
structures), the Hi column density distribution function, the 
nature of the equilibrium between "cold" and "warm" phases, 
etc. In recent times, there have been several indications of a 
surprising degree of small-scale structure in the atomic ISM. 
Perhaps the first of these was the observation of spatially vari- 
able H I absorption towards compact radio sources (e.g. Dieter 
et al. [T9761 Davis et al. [T9961 Faison et al. 119951 . In the 
most extreme cases, for example towards 3C138, there is ev- 
idence for changes of as much as At = 0.1 in the Hi opacity 
on transverse scales of 20 AU (Faison et al. 1998 1. The sim- 
plest interpretation of these observations requires large vari- 
ations in the volume density, Anm ~ 10 5 cirr 3 , assuming 
that all of the other relevant variables (specifically, the path- 
length and the temperature) are kept fixed. However, as ar- 
gued by Deshpande (2000 1, realistic ISM structure functions 
can lead to large variations of r with small angular offsets, 
simply from statistical fluctuations in the effective pathlength 
with position. Further evidence for small-scale atomic structure 
stemmed from searches for time variability in the H I absorp- 
tion seen toward pulsars (e.g. Frail et al. 1994). However, the 
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early claims for ubiquitous and significant time variations in r 
have not been confirmed in more recent studies, (e.g. Johnston 
et al. 2003, Stanimirovic et al. 2003). A third line of evidence 
has come from observations of Na I observation toward nearby 
pairs of stars (eg. Watson & Mever ll99rjl Points et al. [2004 
and references therein), suggesting discrete absorption features 
in the local ISM with highly variable properties on scales of 
hundreds of AU. However, Points et al. (|5004i note that Nal 
is not the dominant ion in diffuse clouds and it is thus possible 
that the above structures are not distinct physical entities but 
merely reflect fluctuations in local physical conditions, such as 
the temperature, pressure, electron density, etc. 

2. Observations 

We have recently undertaken a series of extremely sensi- 
tive H I absorption observations towards bright compact extra- 
galactic radio sources near the North Galactic Pole (NGP), 
with the Westerbork Synthesis Radio Telescope (WSRT). The 
initial motivation for these observations was the detection of 
Hi absorption from a Warm Neutral Medium (WNM), even 
for temperatures as high as 10 4 K. Such an experiment was 
prompted by the detection of wide H I absorption by Kanekar et 
al. (|2003 1, with velocity spreads corresponding to an equivalent 
temperature of 3500 K; these observations reached a l<x opti- 
cal depth sensitivity of ~ 6 x 10~ 4 per 1 km/s velocity channel. 
The NGP region was chosen to ensure that the lines-of-sight 
through the Galactic disk might be as short as possible and, 
therefore, relatively simple. Sensitive, high velocity resolution 
absorption spectra were acquired towards 3C286 (14.7 Jy, 12 
hr integration), 3C287 (7 Jy, 9 hr), 4C+32.44 (5 Jy, 1 1 hr), 
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and B2 1325+32 (1.4 Jy, 11 hr), by observing in an in-band 
frequency-switching mode, utilising a 1 MHz throw every 5 
minutes inside a total bandwidth of 2.5 MHz, with a chan- 
nel width of 0.5 km/s. The in-band frequency switching al- 
lowed exceptionally good band-pass calibration while provid- 
ing 100% of observing time on-source. We achieved RMS op- 
tical depth sensitivities At of 1.8x10 4 (for 3C286), 2.9xl0 -4 
(for 3C287) and 4.0xl0~ 4 (for 4C+32.44) at 1 km/s veloc- 
ity resolution. These are the most sensitive Hi 21cm absorp- 
tion measurements of which we are aware. A somewhat lower 
sensitivity {At ~ 0.0028 per 1 km/s) was obtained towards 
B2 1325+32, due to its lower flux density. Finally, we also used 
the simultaneously acquired total power spectra of the WSRT 
single dishes to derive sensitive H I emission profiles (with ef- 
fective integration times of more than 100 hr) towards all four 
sources. The flux scale is based on the SEFD of the WSRT 
dishes at 1400 MHz of 300+10 Jy/Beam. We emphasise that 
only a constant term has been subtracted from each spectrum 
(both absorption and emission) and not even a first order poly- 
nomial, let alone a higher order spectral baseline. 

3. Results 

Since the primary motivation for the experiment was a search 
for broad, shallow absorption along "simple" high latitude lines 
of sight, we were surprised to find instead that multiple, narrow 
(FWHM ~ 2-3 km/s) absorption features were detected at dis- 
crete line-of-sight velocities, with peak opacities in the range 
0.1-2%. A representative cool phase spin temperature of 100 K 
(corresponding to a thermal FWHM of 2.2 km/s) yields ex- 
tremely low H I column densities, 0.4-8 xlO 18 crrT 2 , for these 
components. The detection of discrete absorption features was 
particularly surprising since the peak emission brightness seen 
in these directions (with the 35'total power beam) is only be- 
tween 2 and 5 K. The absorption and emission spectra towards 
the four sources are compared in the top and bottom panels 
of Figure^ It is striking that moderately nearby lines-of-sight 
(separation < 105' between 3C286 and 4C+32.44) show es- 
sentially uncorrected absorption spectra, and only weakly cor- 
related emission features. Only in the case of 4C+32.44 and 
B2 1325+32, which have an angular separation of only 15', is 
there a hint of similar velocity components having been de- 
tected. In all cases, the emission spectra show wide wings ex- 
tending out to +50 km/s, while the absorption is restricted to 
the velocity range to -30 km/s. The very broadest emission 
in the spectra of Fig.Qb) is actually due to so-called stray radi- 
ation, H I emission attenuated by 30 to 40 dB, but still detected 
in the far side-lobes of the telescope response when the bright 
Galactic plane happens to be above the local horizon. 

The same WSRT observations that provided the absorption 
spectra also allow an interferometric imaging search for emis- 
sion counterparts. Cubes of H I line emission were produced at 
a range of angular resolutions (15, 30, 60 and 120"). Obtaining 
sufficient brightness sensitivity to achieve detections in emis- 
sion typically required smoothing to 60". At those velocities 
where total power emission exceeding about 1 K is seen in 
Fig-H compact emission clumps of 2-3 K brightness are de- 
tected at apparently random locations in the field, superposed 
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Fig. 1. H I absorption (top) and total power emission (bottom) 
spectra after Hanning smoothing to 1 km/s velocity resolution. 

on the poorly-sampled (by the interferometer) diffuse back- 
ground emission. One such emission clump, immediately adja- 
cent to the 3C286 line-of-sight is shown in Fig.|2ja). The intrin- 
sic angular size of these clumps appears to be about 30", while 
their FWHM line-widths are 1-2 km/s, corresponding to tem- 
peratures of ~ 20-80 K. A representative peak column density 
for the clumps in the 3C286 field is N H i ~ 5xl0 18 ctrT 2 . At an 
assumed distance of, say, 100 pc, the clumps would have a size 
of ~ 3000 AU, with a central volume density of ~ 100 cm -3 . 

An even more interesting emission structure is detected in 
the 4C+32.44 field. The line-of-sight toward this background 
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Fig. 2. (a) Top panel: H I emission clump adjacent to 3C286 in 
map (left) and spectrum (right), (b) Bottom panel: H I emission 
shell toward 4C+32.44 in map (left) and spectrum (right). 



source appears to intersect a 15'diameter shell of Hi emis- 
sion, as shown in Fig|2jb). Image quality is limited in Fig|2jb) 
by a non-ideal configuration and no primary beam correction 
has been applied. Although it may be a chance superposi- 
tion, this apparent shell contains the GOIII star HD 116856 at 
(a, (5)2000 = (13:25:55.835,+31:51:40.629). The measured par- 
allax of this star places it at 105 + 11 pc, where the shell 
would have a diameter of ~ 0.45 pc. The stellar proper mo- 
tion (Aa, A6) = (+14.73,-43.40) mas/yr is directed to the SE. 
PeakNm ~ 10 19 cm -2 is seen in this structure, with FWHM 
line-widths of 2-3 km/s. It seems possible that this structure is 
associated with the termination shock of a stellar wind. 



4. Cloud models 

Various approaches are used in the literature to constrain phys- 
ical conditions from an analysis of H I emission and absorption 
spectra (eg. Kanekar et al. 120031 Heiles & Troland 2003). An 
exhaustive comparison of the various methods is beyond the 
scope of the current paper, but will be addressed in Braun & 



Kanekar (2005). We simply note from the outset that while it 
is always possible to model both H I absorption and emission 
profiles as arising from the sum of multiple Gaussians, it is far 
from clear whether a plausible combination of physical con- 
ditions might exist along any real line-of-sight that would pro- 
duce such an artificially simplistic observable (ie. very long, yet 
iso-thermal, and non-turbulent path-lengths). In fact, the most 
sensitive observations of both emission and absorption lines 
(cf. Fig. ^ show profiles that appear to be semi-Lorentzian, 
with narrow line-cores that merge smoothly into broader wings. 
As more sensitive data become available, an increasing number 
of broader Gaussians becomes necessary to fit such profiles. 



Motivated by the non-Gaussian line profiles we have ex- 
plored some simple, spherically symmetric, isobaric two-phase 
cloud models of the form: 



«hO) = n c exp[-(r / sf] 



«hO) = n exp[-(r/sf 2 ] 



forT < 4000 K and 



forT > 4000 K 



(1) 



(2) 



where we further relate volume density to temperature using a 
constant thermal pressure P/ks = n^T = 1500 cm" 3 K. The 
gas temperature was allowed to vary between T, m >, = 20 K and 
Tma.x = 15000 K, yielding an assumed thermal velocity disper- 
sion of cr 2 = 0.0086T. The predicted H I absorption and total 
power emission spectra were calculated for a "cloud" placed 
at the central velocity of each observed feature in an attempt 
to simultaneously reproduce the observed spectra shown in 
Fig-H The most important free parameters in this process were 
the cloud scale-length, s, the cloud distance, d (which most 
strongly influences the predicted total power emission) and an 
impact parameter, b. This last parameter was used to allow for 
the likely circumstance that each spherical model cloud may 
not be penetrated exactly on-axis by the background absorber. 
The two power-law indices of the scaled radius, a\ and a-i de- 
termine the characteristic line shapes in the cold core and warm 
halo respectively. Although these were, in principle, also free 
parameters, it was found that only minor variations from "stan- 
dard" values of about ai=l/4 and a2=l/8 were needed. Our 
arbitrary choice of a transition temperature of 4000 K to sepa- 
rate the cool and warm phases is also non-critical. Comparable 
fits are possible with other choices of the transition tempera- 
ture within a broad range. The broad-band (~ 90 km/s FWHM) 
stray radiation contribution to the emission spectra is approxi- 
mated by a simple Gaussian, rather than a physical component. 

An illustration of the simultaneous emission and absorption 
line fitting is shown in the top panel of Fig. [3] for the 3C286 
line-of-sight. A formal least squares fit has not been carried out 
here, but merely a ^-by-eye to illustrate the possibilities of this 
approach. The bottom panel of the figure shows the entire set of 
density and temperature profiles as well as the fractional distri- 
bution functions for the modelled clouds along all four lines-of- 
sight. This illustrates the basic similarities of the modelled fea- 
tures, although the apparent scale-lengths do vary by about an 
order of magnitude. The apparent distance used in reproducing 
the spectra was about 10 pc, and the typical half-density radius 
was 100 AU, although these were not very well-constrained 
given the very low angular resolution of our total power data 
(35 arcmin). These modelling results can be easily scaled to 
other assumed thermal pressures by a linear scaling of «h to- 
gether with an inverse linear scaling of both s and d. For exam- 
ple, with a typical thermal pressure of only 150 cirT 3 K, com- 
parable fits would be obtained at apparent distances of 100 pc 
and typical half-density radii of 1000 AU. Variations from the 
idealised spherical cloud symmetry would have a similar im- 
pact on apparent distances and sizes. The preserved quantity in 
these scalings is the distribution of column density over tem- 
perature (or more precisely line-width). The fractional column 
density at a particular density and temperature with respect to 
the total line-of-sight column density is plotted in the lower 
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Fig. 3. Overlay of observed, modelled and residual spectra for 
the 3C286 1-o-s (top). (The broad stray radiation contribution to 
the emission profile is plotted separately.) Density and temper- 
ature profiles (left) and fractional distribution functions (right) 
of all modelled "clouds" (bottom). 

right panel of Fig. [3] Column densities of the cool cores are typ- 
ically only a few percent of those in the warm envelopes. The 
total column densities of the modelled clouds vary between 5- 
35xl0 18 cm" 2 . While it should be emphasised that the current 
attempts are fairly simplistic, it is interesting that we obtain 
scale lengths quite similar to those of the structures seen di- 
rectly in emission in Fig|2 We also stress that our cloud model 
"temperatures" may not be wholly thermal in nature, but are 
likely to include a non-thermal, turbulent contribution. 

5. Discussion 

The narrow absorption features of Fig. ^ have the lowest 
column densities that have ever been detected in the CNM, 



two orders of magnitude lower than the mean CNM column 
density Nm = 2.7 x 10 20 cirT 2 of the "classical" McKee- 
Ostriker J1977i model (see also Stanimirovic & Heiles l2005> . 
Our current results thus suggest that even the most diffuse re- 
gions of the Galaxy are populated by hitherto undiscovered, 
tiny distinct structures of very high density and temperature 
contrast. This implies substantial injection of fluctuation power 
on very small scales, whose source is presently unclear. The 
detection of a narrow shell-like structure with small transverse 
separation from a G0III star in one of our fields suggests that 
the stellar winds of intermediate mass stars might play a role 
in the formation of these tiny clumps, whenever such stars 
find themselves within a diffuse atomic structure. More dis- 
tributed sources of energy injection might also lead to ther- 
mal condensation in a turbulent flow as discussed by Audit 
& Hennebelle (2005 1. These authors also demonstrate that lo- 
cal pressure equilibrium is still approximately preserved in the 
vicinity of condensations, despite a wide range of associated 
gas temperatures, consistent with the isobaric two-phase cloud 
models we develop in The density and temperature distri- 
bution functions of our model clouds, Fig.|5]bottom right, also 
compare very favourably with those of Audit & Hennebelle 
(their Fig. 8) in the highly turbulent regime. 

We emphasise that the physical origin of these tiny Hi 
structures and their general importance in the ISM is still un- 
known. Their very low column densities should cause them to 
evaporate rapidly (e.g. McKee & Cowie [T977t . Perhaps the first 
issue that needs to be clarified is whether they are truly ubiq- 
uitous in the ISM. We plan to address this with high sensitiv- 
ity WSRT absorption spectra towards a larger sample of bright 
high-latitude sources. Sensitive total power emission spectra 
with both higher angular resolution and a lower stray-radiation 
contribution using the Green Bank Telescope or Arecibo would 
also be invaluable. Clearly, more work needs to be done to fully 
characterise the nature of sub-structure in the diffuse ISM. 
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